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1. Introduction. Rods are parts of various devices (see [I]- [5]). If it 

is impossible to observe the rod directly, the only source of information 
about possible defects of its fastening can be the natural frequencies of 
O ■ its flexural vibrations. The question arises whether one would be able to 

^ ; detect damage in rod fastening by the natural frequencies of its flexural 

■ vibrations. This paper gives and substantiates a positive answer to this 

CN ■ question. 

The problem in question belongs to the class of inverse problems 
(2Lh'; and is a completely natural problem of identification of the boundary 

\ conditions. 

Closely related formulations of the problem were proposed in P [7]. 
Contrary to this, in this paper it is not the form of the domain or size 
of an object which are sought for but the nature of fastening. The 
^ , problem of determining a boundary condition has been considered in 

^ ! [8]. However, as data for finding the boundary conditions, we take not 

oo ■ a set of natural frequencies, but not condensation and inversion (as in 



[H])._ _ 

O ; Similarly formulated problems also occur in the spectral theory of 

oo , 

O ■ differential operators, where it is required to determine the coefficients 

of a differential equation and the boundary conditions using a set of 
X ; eigenvalues (for more details, see [9]-[l4]). However, as data for finding 

5^ \ the boundary conditions, we take one spectrum but not several spectra 

or other additional spectral data (for example, the spectral function, the 
Weyl function or the so-called weighting numbers) that were used in [9]- 
|13| . Moreover, the principal aim there was to determine the coefficients 
in the equation and not in boundary conditions. The aim of this paper 
is to determine the boundary conditions of the eigenvalue problem from 
its spectrum in the case of a known differential equation. 

The problem of determining a boundary condition using a finite set 
of eigenvalues has been considered previously in [I5]-[I7]. In contrast 
to papers [I5]-[I6], we think it is necessary to determine not the type 
of fastening of plates, but the type of fastening of a rod. The problem 
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of finding of fastening of one of the end of a rod is considered in [TT] . 
In contrast to papers we think it is necessary to determine not the 
type of fastening of one of the ends of a rod, but the type of fastening 
of both ends of a rod. 

2. Formulation of the direct problem. The problem of flexural 
oscillations of rod is reduced to following boundary problem [1| [2] 

[au")" = pFu, (1) 

[(a - q w] = 0, [aw" - C2 = 0, (2) 

[{a u")' + c+ u'] = 0, [a u" + c+ u'] = 0, (3) 

where a is the flexural rigidity, / is the length, p is the density and 
F is the cross-section area of the rod, cf {0 < cf < cx)) are rigidity 
coefficients of springs. 

Let a, p and F be constant. By denote pF ixP' / a. 

For vibrations, we write 

u{x^t) = y{x) cosiot 

and hence obtain the following eigenvalue problem for y{x): 

= s\ (4) 

Ui{y) = -ai y{0) + a, y"\0) = 0, U2{y) = -a^ y'{0) + as y"{0) = 0, 

(5) 

U'i{y) = as y{l) + as y"'{l) = 0, U4{y) = ae y'il) + ay y"{l) = 0, (6) 

where = /oFcj^/a, a^ > = 1, 2, . . . , 8). 

The coefficients a^ characterize conditions for fastening the rod (rigid 
clamping, free support, free edge, fioating fixing, elastic fixing). 

We note that the functions 

yi (x, s) = (cos ^/s X + cosh ^/s a:)/2, 

?/2(x,s) = (sinv^x + sinhv^x)/(2v^), 

y2,{x , s) = {— cos ^/s X + cosh, ^/s x) / {2s) , ^ ^ 

yA{x, s) = (— sin ^/s x + sinh ^/s x) / (2v^), 



are linearly independent solutions of the equation (jll) which satisfy the 
conditions 
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The natural frequencies cui (si) are the corresponding positive eigen- 
values of problem (j4])-(j6]) (see P, E]). The non-zero eigenvalues of prob- 
lem (j4j)-(j6]) are the roots of the determinant 

Ui{yi) U,{y2) Ui{y,) Ui{y,) 

U2{yi) ^2(2/2) U2{y3) U2{y,) 

Usiyi) ^3(2/2) Usiys) ^3(2/4) 

Uiiyi) U,{y2) U,{ys) U,iy,) 

We shall now formulate the direct eigenvalue problem (j4l)-(j6]): it is 
required to find the unknown natural frequencies of the oscillations of 
the rod from = 1, 2, ... 8). 

Thus, finding the natural frequencies cjj = 1, 2, . . . ) is equivalent 
to founding of the roots Sj = 1, 2, . . . ) of A(s). 

Thus, knowing (z = 1, 2, . . . 8) it is possible to find Sj = 1, 2, . . . ) 
by standard methods [SIE]. The solution to this direct problem presents 
no difficulties. The question arises whether one would be able to do 
the reverse and find {i = 1,2,... 8) knowing = 1, 2, . . . ). In 
a broader sense it may be stated as follows. Is it possible to derive 
unknown boundary conditions with a knowledge of s^? The answer to 
this question is given in the next section. 

3. Formulation of the inverse problem. We shall denote the ma- 
trix, consisting of the coefficients (j = 1, 2, . . . , 8) of the forms Ui{y) 
{i = 1, 2, 3, 4) by A and its minors by Mij^m {hj, m, /c = 1, 2, . . . , 8): 

ai 04 
0020300000 
0000 as 00 ag 
ae a? 

In terms of eigenvalue problem (j4l)-(j6]), the inverse problem con- 
structed above should be formulated as follows: the coefficients aj {j = 
1, 2, ... 8) of the forms U^{y) {i = 1, 2, 3, 4) of problem (^-(^ are un- 
known, the rank of the matrix A to make up these coefficients is equal to 
four, the eigenvalues Sk of problem (j4j)-([6]) are known and it is required 
to find the matrix A (the class of linear equivalent matrixes). 

4. The duahty of the solution. Together with problem (|4l)-(j6]), 
let us consider the following eigenvalue problem 

y^'^ = s\ (9) 



A(s) = 



ijkm 



= ± ttj ttj at CLj 
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Ui{y) = -aiy{Q) + a^y"\Q) = 0, U^iy) = -a2 ^/'(O) + ^/"(O) = 0, 

(10) 

UM = «5 y{l) + as y"{l) = 0, U,{y) = ae y\l) + y"{l) = 0. (11) 

^We denote the matrix composed of the coefficients of the forms 
Uj{y) by A and its minors by Mijkm- 



A = 



ai 













a4 



















CLq Oj'j 






OS 





M 



ijkm 



± tti ttj ttk a 



m- 



Let A* be matrix 



a5 as 

06 ar 

ol 04 

02 03 



Theorem 1 (on the duality of the solution of the inverse problem). 
Suppose the following conditions are satisfied 



rank A = rank A = A. 



(12) 



If the eigenvalues {si] of problem (I4l)-(j6]) and the eigenvalues {s^} of 
problem (Ml-dTTl) are identical, with account taken for their multiplic- 
ities, the classes of linearly equivalent matrixes A and ^4 or ^4 and A* 
are also identical. 

Proof. The eigenvalues of problem (j4j)-(j6]) are the roots of the deter- 
minant 



A(s) = -Mi256 [/- + (M2457 + Mi368) [/"(s)] + M3478 [s^ f~{s)] 
+ (Mi278 + M3456 - M2457 - Mi368) [/+(«)] + M1357 [z{s)/s] - M2468 [s z{s)] 
+ (M2456 + M1268) [g"{s)/yM - (^3468 + M2478) [V^ 9~ {s)] + 
+ (Mi356 + M1257) [g+{s)/V?] - (M3457 + Mi378) iV^g^is)], 

(13) 

where 

/^(s) = (1 ± cos ^/s cosh v^)/2, z{s) = (sin sinh ^/s)/2, 
5'^(s) = (— sin y/s cosh \^ ± cos y/s sinh y/s)/2. 
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Note that the functions in square brackets are Hnearly independent. 
Since A(s) ^ 0, A(s) ^ are entire functions in s of order 1/2, it 
follows from Hadamard's factorization theorem (see [18]) that determi- 
nants A(s) and A(s) are connected by the relation 

A(s) = CA(s), 

where /c is a certain integer and C is a certain non-zero constant. From 
this and we obtain the equalities 





= K Mi256, 


(14) 


Mi357 = 


= i^Mi357, 


(15) 


^2468 = 


-i^M2468, 


(16) 


M3478 = 


= i^M3478, 


(17) 


M1257 + Mi356 = 


^-(71^356 + ^1257), 


(18) 


Ml 268 + M2456 = 


A-(Mi268 + M2456), 


(19) 


Mi378 + M3457 = 


^-(71^378 + ^3457), 


(20) 


M2478 + M3468 = 


A-(M2478 + M3468), 


(21) 


Mi278 + M3456 = 


A-(Mi278 + M3456), 


(22) 


Mi368 + M2457 = 


K (Ml368 + M2457). 


(23) 



Hence, 

Cll = ttl, (22 = fl2, 0^3 ~ 0^3, CL4 = CL^^ 

CL^ = OL^^ ClQ = ClQj CL-j = OL-j^ CL^ = (28, 

or 

Gl = 05, CL2 — CLq^ O3 — O7, tt4 = fl8, 

(25 = fll, flg ~ ^4, 0^7 ~ 0^3, (28 ~ ^4- 

These equations can be proved by standard method of algebra. Let 
us consider only one case M3478 7^ for example. If M3478 7^ 0, 03 7^ 0, 
04 7^ 0, 07 7^ 0, 08 7^ 0, 03 7^ 0, (i4 7^ 0, (?7 7^ 0, (is 7^ and so matrixes 
A and A have the forms 

ai 1 
02 100 000 
000 as 001 
ae 1 
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ai 1 

02 100 000 

000 05 001 

og 1 



From this and from ( ITTl ) we obtain K = 1. 
It follows from ([I5]), ([20]) that 



ai as = ai as, ai + as = ai + a^. 



Therefore by Vieta's theorem, we have 



tti — ai. 



as = as 



or 



ai — as. 



as = ai. 



(24) 



Similarly, equations (|T6|l, (EI]) imply 



a2 — a2. 



ae — ae. 



or 



0^2 — 0^6? 0^6 — 0,2- 



(25) 



These equations are equivalent to four cases: 



Case 


1 


ai 


= ai, 


as 


= as. 


a2 


= 02, 


ae 


= 0(3 


Case 


2 


ai 


= 05, 


as 


= Ol, 


a2 


= ^6, 


ae 


= 02 


Case 


3 


ai 


= oi, 


as 


= 05, 


a2 


= ag. 


ae 


= 02 


Case 


4 


ai 


= 05, 


as 


= ai, 


a2 


= 0^2, 


ae 


= Oq 



Cases 3 and 4 are special cases of 1 and 2. Indeed, in case 3 it follows 



Hence ai — as = or a2 — ae = 0. If ai — as = 0, then case 3 is special 
case 2: oi = ai = as = as- If a2 — ae = 0, then case 3 is special case 1: 
ae = a2 = ae = ae- 

In the same way, we can show that the case 4 is special case of cases 
1 and 2. We have proved that in the case M3478 7^ it is possible only 
two cases: 1 and 2. In a similar manner we can prove the analogous 
results for other cases. This contradiction proves the theorem. 

5. Exact solution. It has been shown above that the problem 
of finding the matrix A (in the sense that class of linearly equivalent 
matrixes of ^4) from the natural frequencies of flexural oscillations of a 
rod has a duality solution. The next question is how this solution can 
be constructed. 

This section deals with solving this problem and constructing exact 
solution by the first 9 natural frequencies Ui. 

Suppose si, S2, . . . , sg are the values corresponding to the first nine 
natural frequencies cui. We substitute the values si, S2, . . . , sg into ( 1131 ) 



from {221) or (ESI) ^^at 




(26) 
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and obtain a system of nine homogeneous algebraic equations 



in the ten unknowns 

Xi = -Mi256, X2 = M2457 + ^1368, = M3478, 

X4 = Mi278 + ^3456 " ^2457 " ^1368, 2:5 = M1357, 

Xq = -M2468, = ^2456 + M1268, ^8 = -(M3468 + M2478), 

Xg = Mi356 + M1257, Xio = -(M3457 + M1378). 



The resulting set of equations ( 1271 ) has an infinite number of solutions. 
If the resulting set has a rank of 9, the unknown Xi can be found in 
accurate to a coefficient. The unknown matrix A is found from ( 1311 ) by 
direct calculations as in proof of theorem 1. 

These reasons prove 

Theorem 2 (on the duality of the solution of the inverse problem). If 
the matrix of system ( 1271 ) has a rank of 9, the solution of the inverse 
problem of the reconstruction boundary conditions (ITOl). ( ITTI ) is duality. 

Remark. Theorem 2 is stronger than theorem 1. Theorem 2 use only 
9 natural frequencies for the reconstruction of boundary conditions and 
not all natural frequencies as in theorem 1. But theorem 1 is proves 
duality of the solution in the common case (if the set ( 1271 ) has not a 
rank of 9). 

Continuity of the solution of the inverse problem with respect to Si is 
proved as in [16]. This shows that small perturbations of eigenvalues Si 
{i = 1,2,3) lead to small perturbations of the boundary conditions. It 
follows from this and theorem 1 that the inverse problem is well posed, 
since its solution exists, is unique and continuous with respect to Sj 



Computer calculations confirm the stability of the solution of the 
inverse problem. The order of error often hardly different from the error 
in the closeness of values of Jk and Si and only in some cases it can be 
deteriorated by four orders of magnitude. So the measurement accuracy 
of instruments to measure natural frequencies must exceed accuracy to 
measure boundary conditions by four orders of magnitude. 



+ ^2 [f'{s^)] + X3 [sj r(s,)] 



(27) 




I = 1,2,. ..,9 




(^ = 1,2,3). 
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3. Numerical results. We use dimensionless variables in the nu- 
merical examples. 

Example 1 (rigid clamping — free support). Suppose 

si = 15.4182057169801, S2 = 49.9648620318002, 
S3 = 104.247696458861, S4 = 178.269729494609, 
S5 = 272.030971305025, sq = 385.531421917553, 
S7 = 518.771081332259, sg = 671.749949549144, 
S9 = 844.468026568208 

correspond to the first 9 natural frequencies uji determined using instru- 
ments for measuring the natural frequencies, then the solution of set 
( 1271 ) has the form 

Xi = -Mi256 ~ • C, X2 = M2457 + Ml368 « • C, 

Xi = M3478 « • C, x^ = Mi278 + M3456 - M2457 - Mnm « • C, 

X5 = M1357 « • C, X6 = -M2468 « • C, 

X7 = -M2456 + M1268 ~ • C, = -(M3468 + M2478) « • C, 

X9 = Mi356 + M1257 ~ C, Xio = -(M3457 + Mi378) « • C 

(29) 

with an accuracy of 10 ^ . 

Substituting 1 for C in (1291), we get M1356 + M1257 = 1. This means 
that Mi356 7^ or M1257 7^ 0. 

Suppose Mi356 7^ 0. Then ai 03 05 7^ 0. So 

1 a4 

02 10000 
0000100 ag 

000001 07 

(in the sense that class of linearly equivalent matrixes of ^4). 

Using Mi256 = 0, M1357 = and M2457 + ^1368 = 0, we have 02 = 0, 
07 = and ag = 0. From this and Mi27g + ^3455 — M2457 — Mi36g = 0, 
we get a4 = 0. 

Thus we have 

1 
00100000 
00001000 
1 
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Suppose now M1257 7^ 0. Then ai 02 05 a? 7^ 0. So 



1 a4 

1 as 

1 as 

as 1 



Using Mi256 = and M2457 + M1368 = 0, we get ae = and 04 = 0. 
If we combine this with equations M1357 = and M1278 + ^3455 — 

M2457 — Mi368 = 

Thus we have 



we get 03 




and 


as 




0. 






1 




















A = 





1 






























1 
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So boundary conditions are or 
^i(2/)=?/(0) = 0, 

Usiy) = y{i) = 0, 



either 



Ui{y) = y{0) = 0, 

Uz[y) = y{l) = 0, 



U2{y) = y"{0) 
U,{y) = y'il) - 

U2{y) = y'iO) 
U,[y) = y"[l) 



--0, 
0, 

0, 
0. 



Note that the numbers Si presented above are the same as the first 
nine exact values corresponding to fastening of rigid clamping and free 
support. This means that the unknown rod fastening inaccessible to 
direct observation has been correctly determined. 

Example 2 (rigid damping — elastic fixing). If 

si = 5.60163863016235, S2 = 22.4984332740862, 
S3 = 61.8604321649037, S4 = 120.984868139371, 
S5 = 199.909638169628, = 298.589053349029, 
S7 = 417.014779762035, sg = 555.183266366176, 
sg = 713.092945010199 

correspond to the first 9 natural frequencies cui determined using instru- 
ments for measuring the natural frequencies, then the solution of set 
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( 1271 ) has the form 



Xi = -Mi256 ~ 0, X2 = M2457 + Mi368 ~ 0, 

= Mm78 ~ 0, X4 = Mi278 + ^3455 - M2457 " ^1368 ~ "C, 

X5 = M1357 « 0, xq = -M2468 ~ 0, 

X7 = -M2456 + M1268 « 0, X8 = -(M3468 + M2478) ~ 0, 
X9 = Mi356 + M1257 ~ 5 C, Xio = -(M3457 + ^1378) « 0. 

(30) 

Suppose C = 1; then M1356 + M1257 7^ 0. This means that M1356 7^ 
or M1257 7^ 0. 

Suppose Mi356 7^ 0. Then ai 05 ae 7^ 0. So 

100 04 000 
02 10000 

0000100 08 

1 ay 

(in the sense that class of linearly equivalent matrixes of A). 

If we combine this with equations M1256 = 0, M1356 + M1257 = 5 C we 
get 02 = 0, 1 = 5 C. Using M1357 = and M2457 + M1368 = 0, we have 
aj = and = 0. From this and M1278 + M3456 — M2457 — M1368 = — C, 
we get —04 = —1/5. 

Thus we have 

1 1/5 
001 0000 
000 1000 
1 

Suppose now M1257 7^ 0. Then ai 02 <^5 cl7 7^ 0. So 

100 a4 0000 

01 a3 00000 
000 0100 a8 
000 OOaelO 

Using Mi256 = and M2457 + M1368 = 0, we get aQ = and 04 = 0. 

If we combine this with equations M1357 = 0, M1356 + M1257 = 5 C, 
Mi278 + M3456 - M2457 - ^1368 = -C and M2457 + ^1368 = we get 
03 = 0, 1 = 5 C, and —ag = —1/5. 
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Thus we have 



1 

0100000 

1 1/5 

1 

So boundary conditions are or 

Ui{y) = -5 y{0) + y"'{0) = 0, U2{y) = y"{0) = 0, 
Usiy) = y{l) = 0, U,{y) = y\l) = 0, 

either 

Ui{y) = y{0) = {), U2iy)=y'i0) = 0, 
UM = ^y{l)+ y"\l) = 0, U,{y) = y"{l) = 0. 

Note that the numbers Si presented above are the same as the first 
nine exact values corresponding to fastening of rigid clamping and elastic 
fixing with the relative stiffness factor of 5. This means that the un- 
known rod fastening inaccessible to direct observation has been correctly 
determined. 

Example 3. (elastic fixing — elastic fixing) If 

si = 0.383848559322840, S2 = 11.9180148367849, 
S3 = 53.4326824121208, S4 = 114.157790468867, 
S5 = 193.836586296759, se = 292.955617117120, 
S7 = 411.667770695782, sg = 550.037492353361, 
S9 = 708.096219400352 

correspond to the first 9 natural frequencies ui determined using instru- 
ments for measuring the natural frequencies with an accuracy of 10^^^, 
then the solution of set ( |27l ) has the form 

XI = -Mu56 « -24 C, X2 = M2457 + Mi368 « "10 C, 

X3 = M3478 « C, X4 = Mi278 + M3456 - M2457 - Mi368 « 9 C, 

X5 = M1357 « 3 C, xq = -M2468 ~ -8 C, 

X7 = -M2456 + M1268 « -32 C, xs = -(M3468 + M2478) « -6 C, 

XQ = Mi356 + M1257 « 18 C, Xio = -(M3457 + Mi378) « 4 C. 

(31) 
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Suppose C = 1; then M3478 = 03040708 7^ 0. So 



01 1 

02 1 

05 1 

06 1 



(in the sense that class of Hnearly equivalent matrixes of ^4). 

From this and (ISTl). we obtain or oi = 1, 02 = 2, 05 = 3, oe = 4 
either oi = 3, 02 = 4, 05 = 1, oe = 2. So boundary conditions are or 



^i(2/) = -?/(0) + ?/'''(0) = 0, 
t/3(?/) = 3 2/(0 + = 0, 

either 

^i(2/) = -3?/(0) + 2/'"(0) = 0, 

^3(2/) = 2/(/) + y"\i) = 0, 



/72(2/) = -2 2/'(0) + 2/"(0) = 0, 
^4(?/) = 4 ?/'(/) + 2/"(/) = 0, 

^2(2/) = -4 2/'(0)+2/"(0) = 0, 
/74(2/) = 2^'(/) + /(/) = 0. 



Note that the numbers Si presented above are the same as the first 
nine exact values corresponding to elastic fixing with the relative stiff- 
ness factors of 1, 2, 3, 4. This means that the unknown rod fastening 
inaccessible to direct observation has been correctly determined. 
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